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WEAK ERROR FOR CONTINUOUS TIME MARKOV CHAINS 
RELATED TO FRACTIONAL IN TIME P(I)DES 


M. KELBERT, V. KONAKOV, AND S. MENOZZI 


Abstract. We provide sharp error bounds for the difference between the tran¬ 
sition densities of some multidimensional Continuous Time Markov Chains 
(CTMC) and the fundamental solutions of some fractional in time Partial 
(Integro) Differential Equations (P(I)DEs). Namely, we consider equations 
involving a time fractional derivative of Caputo type and a spatial operator 
corresponding to the generator of a non degenerate Brownian or stable driven 
Stochastic Differential Equation (SDE). 


1. Introduction 

are interested in the probabilistic approximation of P(I)DEs of the following 

jd^u(t,x) = Lu(t,x), (t, x) £ R)j_ x R d , 

|u(0,a;) = f(x), x £ 

df,/3 £ (0,1), stands for the Caputo-Dzherbashyan derivative and L is the 
generator of a Brownian or stable driven non degenerate SDE (see equations (12. 101) , 
(12.111) below for the respective definitions of the Caputo-Dzherbashyan derivative 
and the spatial operators considered). Equations of the previous type appear in 
many applicative fields from natural sciences to finance, see e.g. Meerschaert et al., 
[MSI 2) . [MS 131 and references therein. 

Under suitable assumptions, the solution of m can be represented as u(t, x) := 
E[f(X°’p)] where (A°’ a: ) s >o solves the SDE with generator L and (Zf) t >o is the 

inverse of a stable subordinator of index /3 £ (0,1) independent from X 0,x . This 
therefore extends the “usual” Feynman-Kac representation, corresponding to /3 = 1, 
to the fractional case (3 £ (0,1). 

Many probabilistic numerical approximations of u(t, x) have been considered 
when (3 = 1. We can for instance mention the works of Konakov et al. (see [KMOOj . 
[KM02] , IKM09I in the non degenerate diffusive case or [KM101 for SDEs driven by 
symmetric stable processes) that investigate the Euler scheme or, more generally, 
the Markov Chain approximation of the spatial motions in terms of Edgeworth 
expansions or Local Limit Theorems (LLTs). We can also refer to the works of Bally 
and Talay for the Euler scheme of some hypoelliptic diffusions IS , [BT96bj or 
[KMM10] for associated LLTs. 


We 

type: 

( 1 . 1 ) 

where 
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In the current (strictly) fractional framework two additional difficulties appear. 
Firstly, it is known that the fundamental solutions to EH exhibit, additionally 
to the usual time singularity in short time, a diagonal spatial singularity, see e.g. 
Eidelman and Kochubei KKO 1 for the current case or Kochubei |Kocl4) for ex¬ 
tensions to higher order fractional derivatives /3 G (1,2). Secondly, the inverse of 
the subordinator leads to consider random integration times that might be either 
very small or very long. The analysis of the discretization error in the previously 
mentioned works was always performed for a fixed final time horizon T and the 
constants controlling the error estimates depend explosively on T for small and 
large times. We must therefore carefully control these explosions. In short time 
we must handle the spatial singularity of the fractional heat kernel whereas in long 
time, the explosion is compensated by the fast decay of the density of the inverse 
subordinator. Let us mention that those difficulties may deteriorate the “usual” 
convergence rate for the weak error even for the Euler approximation. We estab¬ 
lish error bounds for the Euler scheme and the Markov Chain approximation of 
a diffusive SDE, Theorem O and for the Euler scheme of a stable driven SDE, 
Theorem l3.2l We emphasize that Kolokoltsov [Kol09] also considered probabilistic 
schemes to approximate equations of type EH in a more general setting, namely 
considering a fractional like time derivative that could depend on space as well. He 
establishes convergence of the schemes towards the expected solution but since the 
approach relies on semigroup techniques, no convergence rates are provided. 

The article is organized as follows. We first recall in Section [2] which are the 
probabilistic tools needed to give representations and approximations of the solution 
to equation EH). We then state our main results in Section [3] The proofs are 
given in Section [|] which is the technical core of the work. Some perspectives are 
considered in Section [5j Technical results concerning the parametrix expansions, 
which are needed to establish the short and long time behaviour of the involved 
spatial densities, are collected in Appendix [A] As a by-product of our analysis, 
we obtain two-sided heat kernel estimates for the fundamental solution of EH in 
Appendix Q3] 


2. Probabilistic objects associated with Cauchy Problems with Time 

Fractional Derivatives. 

We first recall how a probabilistic representation for the solution of EH can 
be derived considering the special case of L being associated with a d-dimensional 
symmetric stable process of index a G (0,2], thus including the Brownian motion. 
We write for </> G CQ(R d ,R) (space of twice continuously differentiable functions 
with compact support), Lcj )( x ) = \ A ( j ){ x ) for a = 2 whereas for a G (0, 2): 

L ( f >{ x ) = / {( j){x + y ) - c /)( x ) - V</>(x) • fi { dy ) 

JR d 

(2.1) = [ {(/}(x + \ y \ y ) - ( j )( x ) - V0(x) • \ y \ yl \ y \< i )}\ yr (1+a) d \ y \ p ( dy ), 

JR+xS ^- 1 

where v is the spherical part of /j and is assumed to be non degenerate, i.e. there 
exists 


[ i<£,d>rM^)<Aier. 


( 2 . 2 ) 


A > 1, Vf G R d , A -1 |£| a < 


s d - 1 
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Observe that if (S“) u >o is a stable process with generator L then 


Vt > 0, G R d , ■= E[exp(*(^, Sf))] := exp ^-tC a ^ |(£, rf)\°‘v{dri) S j , 

where the explicit value of C a can be found in Sato |Sat05j . 

Let us now discuss heuristically how some suitable scaling limits of Continuous 
Time Random Walks (CTRWs) actually give a probabilistic interpretation of the 
solution to (ED- These objects thus provide a natural approximation scheme to 
m- Basically, CTRWs are random walks that wait a certain amount of time 
between their jumps. In the i.i.d. case, the mechanism can be described as follows: 
let (Tj, be a sequence of M + x Revalued i.i.d pairs of random variables 

defined on some probability space (f2, J 7 , P). For all t > 0 one defines: 

Nt m 

r t := y ^Yj, N t := max{m G N : ^ T, ; < t}. 

2 =0 2=1 

If the (T,)i 6 N* are exponentially distributed and independent of the (Y/)j e N* the 
process (r t ) t >o is easily shown to be Markovian. This property clearly fails in the 
general case. Of particular interest is the case when the (Ti)i 6 N* and the 
are independent and respectively belong to the domain of attraction of a /3-stable 
and a-stable law with /3 G (0,1) and a G (0, 2]. This property indeed yields that : 

[ nu\ [nu\ 

(2.3) Ti) u > 0 =► (S * + ) u > 0 , and (n“ 1/a £ ^)«>o =► (S£)„>o, 

2=1 2=1 


where |_‘J stands for the integer part and => indicates the usual convergence in 
law for processes. We also used the convention that Ei=i = 0- I n (EO), 5 /3,+ 

is 

a /3-stable subordinator, i.e. a Levy process with positive jumps and its Laplace 
transform writes i/j s p, + (\) = E[exp(— \S & ,+ )] = exp(— uX 13 ), A > 0. On the other 
hand, we will assume that S a is a symmetric Revalued stable process with gener¬ 
ator as in ED- 

Now, the rescaled process associated with the number of jumps (N t )t>o also has 
a limit, namely (n _/3 /V J lt )t>o ~> {%t )t>o, where Zf := inf{s > 0 : S^ ,+ > t} which 
is the inverse process of a /3-stable subordinator. Since ,+ is increasing in s, zf 
also corresponds to the first passage time of /E ,+ above the level t. Thus, one 
formally has: 


(n-^“T nt ) t > 0 = (V)- 1/a 


t> 0 


/ n^(n ^ N nt ) N 

= (n^ 3 )- 1 ^ E Yi 


t> 0 


te) 


t> 0 


From the independence of S a and the limit random variable has, for any 

z t 

t > 0, an explicit density that writes: 

/• + 00 

(2.4) 


r~\-oo 

q(t,x)= ps a (u, x)p Z ft (t, u)dh 
Jo 


where ps^(u,x) = f Rd exp(— i(x, \))<ps% (A)dA and p z p(t,u) stands for the 
density of Zf at point u. Observe that since P[Zf < u\ = P[,S'f’ + > t] we have 
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the following relation between the density of the inverse ZP and the density of the 
subordinator 5 ,/3,+ : 

(2.5) p z «{t,u) = <9 U P [Zf < u] = d u (l -P[Sf’+ < t}) = -d u [ p s p,+ (u,y)dy. 

Jo 

We also have the following important equation for this density. 

Proposition 2.1 (PDE associated with pze)- For fixed t > 0, the density p z p(t,u) 
satisfies in the distributional sense the equation: 


( 2 . 6 ) 


3>f + d u )p Z f>(t,u) = 


t~P 


-6{i 


F(1 — (3) 

where S stands for the Dirac mass at 0 and lD>f denotes the Riemann-Liouville 
derivative. Namely, for a function h and /3 £ (0,1): 


®?h(t) = 


1 


/ h(t — u)u Pdu 

Jo 


r(i - (3) dt U o 

Proof. Let us differentiate under the integral in (12.51) to get for u > 0: 

( 2 - 7 ) 

Pzf>(t,u) = - d u p S (>,+ (u 1 y)dy = — / L*’Pp sf >,+ (u,y)dy = -L*’P p sf >,+ (u,y)dy , 
Jo Jo Jo 

using the Kolmogorov forward equation for the density of the subordinator and 
denoting by L*’P the adjoint of its generator. Precisely, 

1 r+°° 

Recall now that for a smooth (say C 1 ) function <p : R + 
on R _ , for all t > 0, 


ds 

s 1+ P' 

that is extended by 0 




l 


F(-/3) 

1 

rH) 

i 


l Mt - s > - - w) 


tp(t) - ip(0) 


PtP Jo * {t S)S ft fttP / 


F(-/3)(-/3) 

Thus we get the result for u > 0 from m and m differentiating w.r.t. u. 
The proof can be achieved observing that pgp,+ (u,y)dy —> 5(dy) which implies 

u —>-0 

lim u _ ) . 0 + p z p{t,u) = Df$(f) = rhzim : #(i) = It>o- This observation implies the 
form of the singularity at the origin in Equation (12.61) . □ 

We refer to IPod99] for additional general references concerning fractional cal¬ 
culus. From (12.41) and (12.61) we therefore derive that for all t > 0: 

r+oo £-/3 

= J PS<* (u, x){-d u p z p ( t , u) + ^ 6{u)}du 

r+oo 

= / d u psi(u,x)p Z fi(t,u)du+ —-xv<5(a0 


F(l-/3) 


t~P 


T(1 — P) 


$( x )> 


( 2 . 8 ) 


Lq(t, x ) + 
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where we used that d u ps a (u, x ) = Lps<* (u, x) (the spatial density satisfies the 
Kolmogorov backward equation) and the fact that 5“ —> 0. We have also denoted 

u—>■ 0 

by 6 , with a slight abuse of notation, the Dirac mass at 0 in M d . 

In other words, the transition p.d.f. q of (5“ 3 )t>o can be viewed as the funda- 

mental solution of (12.81) . Hence, for a given continuous initial data /, the associated 
Cauchy problem solved by u(t, x) = E [f(x + 5“^)] writes: 


D fu(t, x) = Lu(t, x) + r(i-p) f( x )i x ) e +* 
(t(0,cc) = f(x), x £ R d . 


Alternatively, equation (12.91) can be rewritten in terms of the Caputo-Dzherbashyan 
fractional derivative, denoted df hereafter, recalling that for a smooth function h 
and /? £ (0,1): 


( 2 . 10 ) 


dt h (t) ■■= h'(t - u)u p du , 

D fM*) ^ T(1 _ I3) h ^ = d t h ®’ 1 > °' 


The probabilistic representation of the solution to equation m is thus derived 
from dUSD applying identity (12.101) to the function h(t) = u(t,x). 

Let us mention that, in the special case L = A 2 A, some explicit expansions of 
the density of \/2XS ]L have been obtained by Beghin and Orsingher |BQ09j . 

The correspondence between the previously described CTRWs and Equation 
(12.81) suggested various extensions. The first one consists in introducing a Markov 
Chain for the spatial motion, i.e. the * th -spatial jump depends on the position Yj_i- 
This procedure can lead to consider at the limit differential operators of the more 
general form 

L\4>{x) = (b(x), V^>(x)) + ~Tiia<j*{x)Dl(j)(x)) or 

L 2 4>(x) = (b(x),V(/>(x)) 

+ / {</>( x + o-(x)y) - (f>(x) - Vct>(x) ■ a(x)yli r7 ( x ) y \< 1 }M(dy), 
jR d 

( 2 . 11 ) 

for p(dy) := \y\~ {1+a) d\y\v(dy) as in EB, in equations m , (tnt - We consider 
homogeneous operators in time for notational simplicity. Let (A s ) s >o solve the 
SDE with generator given in (12.111) . i.e. 


b(X s )ds+ [ a(X s ~)dY s , 

Jo 

where Y is a Brownian motion in the diffusive case and a symmetric a-stable 
with a £ (0,2) process with Levy measure /i otherwise. Then, under suitable 
assumptions, needed to get smoothness of the p.d.f. for the spatial motion (X t ) t> o, 
we can derive similarly to the previous computations that the solution to the Cauchy 
problem d writes: 

(2.13) W(t,x) £ R+ x R d , u(t, x) = E[f(X°J)\. 

Z t 


( 2 . 12 ) 


X t = 
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Observe that, similarly to equation (12.41) . denoting by px p 

Z T 


time T > 0, one has for all (x,y) £ (R d ) 2 : 


the density of X z p at 


(2.14) 



p{u,x,y)p zf i ( T,u)du , 


where p(u,x ,.) stands for the density of X° ,x . Let us also mention that modifying 
the laws of the waiting times including a possible dependence on the previous jump 
time and spatial position, leads to consider at the limit a modified fractional deriv¬ 
ative in time. We refer to Kolokoltsov IKolOOl for details and results concerning the 
extension of the formula (12.1311 and convergence of the associated approximations. 


3. Assumptions and Main Results 

3.1. Assumptions on the Coefficients and Approximation Scheme. We as¬ 
sume that the coefficients in equation (12.121) satisfy the following conditions: 

(S) The coefficients b : K. d —> R. d , a : R d —> R d <g) R d are assumed to be bounded as 
well as their derivatives up to order 6. 

(UE) There exists A > 1 s.t. for all (x,£) £ (R d ) 2 : 

A -1 |£| 2 < (crcr*(a:)£,£) < A|£| 2 . 

Now, we consider for a given time step h > 0, setting for i £ N, U := ih , the 
following Markov Chain approximation of equation (12.121) : 

(3.1) W £ N, X£ +i = X* + b{X*)h + a{Xl)h l / a m+1 , X$ = x, 

for i.i.d. Revalued random variables (r]j)j> i and a £ (0,2]. 

Also, since we are considering the fractional derivative, we can explicitly and 
exactly simulate (5 , ^’ + )igN, see, e.g., Weron and Weron [WR95 i. From the exact 
simulation on the time grid, a natural choice to approximate the inverse process 
zf is to consider the discrete inverse process: 

(3.2) Vi > 0, Z?’ h := inf{Sj := ih : Sf;+ > *}. 

Thus, for a given T > 0, we finally approximate X z p appearing in (12.131) by 

vh 

^Z§,’ h ' 

Let us mention that the smoothness assumption (S) could be weakened in the 
case of the Euler scheme. These smoothness conditions are actually those required 
to apply the results of [KM09I . 

3.2. Diffusive case: This case is associated with spatial motions with generators 
of the form L\ in (12.111) . We then specifically assume that a = 2 in (13.11) and that 
the (r/j)j>i are s.t. their moments coincide with those of the standard Gaussian 
law in up to order 2. We consider two kinds of conditions: 

(Ieui) The i have standard Gaussian densities (so that the above scheme is 

the usual Euler discretization). 

(I m ) For a given integer m > 2 (d + 1), the {r]j)j> i have a density Qm which is 
C 4 and has, together with its derivatives up to order 4, polynomial decay of order 
M > d(2m + 1) + 4. Namely, for all z £ R d and multi-index u, \v\ < 4: 

\D v Q m (z)\ < C m (1 + \z\)~ m . 
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We say that assumption (Ad,eui) (resp. (A D,m)) is in force as soon as (S), (UE), 
(Ieui) (resp. (I m )) hold. We write (A d) whenever (A,d,eui) or (A D,m) holds. 

It is well known that under (A^,) the random variables X t ,X^. have densities for 
all t > 0, i e N* respectively, see e.g. IKMOO ;. This property then transfers to X z p 

and X^ 3 h thanks to the convolution equation (12.1411 and its discrete analogue for 
the approximation (see equation (13.31) below). We will denote by px « and p X h 

Z T zP’ h 

the associated p.d.f. To distinguish precisely the approximations under (A/yEui) 
or (Ap m ) we will specifically denote, when needed, by X^' h ,p y Eu\,h the Euler 

T 

approximation and its density. 

We now have the following convergence result: 


Theorem 3.1 (Error Bounds for the Approximation Schemes). 

- Euler Scheme. There exists c := c(( Ap)) s.t. for a given time step h £ (0,1), 
a deterministic time horizon T > 0 s.t. T > h 1 ^ 2 (not necessarily a multiple of h), 
for every i/j/we have under (Au^ui) (Euler scheme): 


\Px« (x,y) - p x ™,h(x,y)\ 


p ,h 


(3.3) = | [ p z e(T,u)p(u,x,y)du-y2P[ZP’ h = t i ]p Eul ’ h (t i ,x,y)\ 

Jm+ i>x 


^ eh f £'/3,Time(T, X y^j T Space (T, X I/) 


where: 

£/3,Tim e(T, X 7/) : = 

^/3,Space(^: X y ) . 

with 

p 0 (T,x-y) := 


(3.4) p@(T, x — y) := 


( Id< 2 Id>3 

\TP/ 2 \x-y I \x — y\ 2 

Id=l + Id>3 . Id=2 \ 


pp(T,x-y) + 7 j^p 0 (T,x-y) 


\x ~ y | 
Id< 2 


TP/ 2 \x - y l 11 ^^ 

- i\x~y\ 2 


+ x - y) + x - y)}> 

+ exp(c Tm) 


x exp [ —c 


TP 


exp (cT^/( 1+/3 )) 
TP d / 2 


■ exp i — c 


-i 


■~y ? 

TP 


2 V(2-/3) 


The two auxiliary functions p,pp are s.t. their sum serves as an upper bound for 

both p x », p x Em,h. Precisely, there exists c := c(/3) s.t. for all T > 0, ( x,y ) £ 
Z T Z P’ h 

(K d ) 2 , x ^ y: 

{px- p +p x Eni,h)(x,y) <c(pp+pp)(T,x-y). 

Z T z p,h 

We refer to Corollary 14.11 for a proof of this last statement. We also establish in 
Appendix iBl a similar lower bound, proving the estimate is sharp. 
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- Markov Chain. Consider now assumption (A D,m) (general Markov Chain). 
Then, there exists c := c((Aiy m )) and for all e £ (0,1/5), c e := c £ ((An,m)) s.t. 
for all given time step h £ (0,1), T& > /i 1 / 5_e , (x,y) £ (R d ) 2 , x ^ y: 

\px« (x, y) - p x » n „ (x, y)\ < c{h£p^ ime {T, x - y) 

(3.5) T z - 

+C e {ft 1/2 ^fspace,LLT(^ X — y) + ^Space.NoLLT^) X ~V> ^)}}> 
where the contribution fg.Tim e ( T,x — y) due to the time sensitivity, see equations 
m and G3D, is as above and: 


pM 

^f3, Space 


,LLT (T,x-y) := | d ^ * ^~ 3 + ^pJ^jVm-[(d-i)+i d=1 }A T i x - v) 


-fm vmA T i x ^v) 


denoting as well for all l £ N*: 


qi,p{T,x — y) := 


I<2< 2 


TP/ 2 \x — y\ Id = 2 


+ w\ exp<cT,,/2) 


X 1 + 


s- 2/1 

TP/ 2 


-i 


(3.6) q m ,p{T,x - y) := 
Also: 

? M,e 


TP d / 2 




1 + 


exp —c 


k~2/| 2 


+ 1 


TP/ 2 

[h}/ b ~ e ^h l d< 2 

* - 2/1 


exp (cTP/( 1+ P' ) ) { \x — y\\ ^ 2_ ' 3 ^ 


f ”s P ace,NoLLT( T ’ X V,h) _ ^- 2 +I d < 2 ex P (cT 0/2 ) 

—2(m—l)+(d—2)+Id<2 


(/lV6-e)l/2 


Comments on the results. 

- Euler Scheme. Observe that the term pp in (13.41) comes from the contribution of 
small times in fomula (12.141) . This gives the additional (diagonal) spatial singularity 
in the fractional case. The term pp in (13.41) comes from the integration over long 
times in (12.141) . It induces a loss of concentration in space w.r.t to the usual case 
P = 1, though preserving the usual parabolic equilibrium: the spatial contribution 
remains comparable to the square root of time. 

• In the diagonal region, i.e. for points s.t. \x — y\ < kTP/ 2 for some k > 1, if we 
restrict to spatial points that are equivalent to the time contribution for the usual 
parabolic metric, i.e. k _1 T^/ 2 < \x — y\ < kT^/ 2 , the previous error bound reads 
in small time, i.e. for T < 1 : 

ch 

(3.7) \Px zP (x,y) -p x Eui,h(x,y)\ < —g {pp + pp){T,x - y). 

T z t 

This means that we find the upper-bound for the density of the processes involved, 
multiplied by a factor corresponding to the singularity induced by a time derivative 
(or equivalently second order derivatives in space) of the non degenerate diffusive 
heat kernel at time . There is even here a difference w.r.t. to the usual anal¬ 
ysis, see e.g. |KM02| . which involves a singularity of order one in space, which in 
the current bound comes from the contribution f/ 3 ,Space(7 1 , x — y) in (13.31) and is 
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negligible if T < 1 . This is due to the discrete approximation of the inverse of 
the subordinator which yields to consider derivatives in time for the densities of 
the spatial motion (see equation (Id.II) and Lemma Id. II below). If the error bound 
(13.31) is now considered for T sufficiently large, the contribution £/3,Spac e(T,x — y ) 
dominates and we are back to the usual error rate for the considered region. 

Observe as well that the spatial diagonal singularity induced by the fractional 
time derivative yields additional constraints. Namely, for the Euler approach, to 
get a convergence rate of order h e one needs to take: \x — y\ > d > 1. 

• In the off-diagonal region, i.e. \x — y\ > nT ^/ 2 , we do not have spatial singularities 
anymore. Indeed, the error bound also reads as in (13771) . Anyhow, the specificity 
comes from the loss of concentration due to the term pp coming from the long 
time integration in (I2.1dl) and appearing in the two-sided estimates of Appendix [Bl 
We refer to Lemmas 14.81 and 14.111 (see also Corollary 14.11 and Theorem IB. II) for a 
detailed presentation of these facts. 

The previous facts make us feel that the previous bounds are sharp, up to con¬ 
stants. 

- Markov Chains. For the error bound (13.51) associated with the Markov Chain 

approximation, the term £ppn me is similar to the Euler case, since it does not 
depend on the chosen discretization scheme but only on the time sensitivities of 
the spatial density, see again (ED. On the other hand two additional contributions 
appear. The first one, £g* Space LLX enjoys the same convergence rate h 1 ^ 2 as in 
the usual LLT (see Petrov [Pet05] . KMOO ] ) and the terms q m -(d-i)+i d =i,p>Qm,p 
appearing in ED are the Markov chain analogue of the contributions pp and pp 
in Ein. Note first that we keep under the same index m as the one 

appearing for the weak error at a given fixed time, see (14.271) below and | !KM09j . 
[KMOOj . for the contribution q m ,p associated with the large time integration and 
which also yields a loss of concentration. Observe as well that the spatial diagonal 
singularities in small time also deteriorate the concentration in this framework. We 
point out that, for the LLT to apply, the considered time has to be sufficiently 
large, namely u > in (12.141) . The last term ^yspace,N oLLT comes precisely 

from those very small times u < /i 1 / 5-6 for which the limit results do not apply. 

- Continuity of the Estimates when /3 f 1- In this case, the Caputo derivative in 

(12.101) tends to the usual derivative. It is therefore natural to ask whether our results 
match the ones previously obtained in that case (see [KM00I . |KM02j . KM l()j S. Let 
us now emphasize that when (3 | 1 then —> exp(—An), A > 0. Thus, the 

subordinator S' ,9 ’ + , and its inverse Z 13 both tend to the deterministic drift with slope 
1. This means that for a given T > 0,P[Z x G du] -—> Spidu). Hence, the time 

integral in (12.141) disappears at the limit, i.e. px p (x,y) —> p(T,x,y). From (13.31) 

z t /3->i 

and ED, the same phenomenon occurs for the scheme provided T = Kh, K G N* 
and the error associated with the time sensitivity, term £\ in (14.11) yielding the 
contribution f^Time would vanish in the previous error bounds. Since the previous 
Dirac convergence would also kill the additional contributions in small and long 
time, we finally derive that the error bounds obtained in (13.31) and (13.51) are coherent 
with those of the previous works at the limit. 

^We recall that (pp + pp)(T, x — y) serves as an upper bound for (px , )(x,y) (see Corollary 

Z T 

mi . 
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3.3. Strictly Stable case. Additionally to assumptions (A), (UE) we will assume 
that : 

(ND) The spherical part v of the measure p, in (12.111) satisfies the non-degeneracy 
condition of equation (IQ) . Moreover, v has a C 3 density w.r.t. the Lebesgue 
measure of the sphere. 

(B) The drift b = 0 if a < 1. 

This last condition is rather usual to have pointwise bounds on the density of the 
SDE, see [KolOOl . We will consider similarly to the diffusive case the approximation 
E3D for a £ (0,2) restricting ourselves to the Euler scheme case, i.e. the i 

are i.i.d. symmetric stable random variables with common law Y\ (driving process 
at time 1). We say that (As) holds if (S), (UE), (ND), (B) are in force. Recall 
from Kolokoltsov [KolOO] and jKMIO j that the random variables Xt,X^ Eul have 
a density for all t > 0, i £ N* respectively. Thus, similarly to the diffusive case, 
this property transfers to the random variables X z p , X d ’ E ^ 1 from (12.141) . With the 
previous notations, we have the following approximation result. 


Theorem 3.2. Under (As), there exists a constant c s.t. for a given time step h £ 
(0,1) and a fixed time horizon T > h 1 ^, and any ( x , y) £ (R d ) 2 s.t. \x — y\ > h x ! a 
we have the following convergence result: 


(3.8) 


ch{ 

LV \x — y\ 


a. 


x-y\>Tf>/ a 

TP 


I Px» (x, y) - p x h ,eui 

Z T z P,h 

)pp(T,x~y) + —7gP/3(T, x - 


(x,y) I 
-y)}» 


where 

Pp(T,x- y) 
Pp(T,x- y) 



exp(cT^ )){ 2 1 /3| x _ y\d-a^\ x -y\< Tf)/a | iT _ y^d+oJ L \^~v\> Tfi,a \' 

exp l c T^/^-^-P))) l 

u : = - A 1, 

TW« (l + “ 

p x Eni,h)(x,y) < c(pp +pp)(T,x — y), c:=c(/3,{A s )). 
zP- h 


Remark 3.1. We observe a phenomenon which is similar to the diffusive case, i.e. the 
error bound has the expected rate of order h up to additional singularities. In the 
strictly stable case, i.e. a £ (0,2), two contributions give additional singularities 
w.r.t. those already appearing in the sharp density bounds (see Theorem IB. ID . 
Those contributions write T -/3 , \x — y\~ a and come from the time derivative of the 
density in (12.141) . 

On the other hand, since the spatial motion already has heavy tails, we do not 
observe here a loss of concentration phenomenon as we did in the diffusive case. 
Eventually, we still have diagonal spatial singularities, which are in the strictly 
stable case more direct to formulate. They appear precisely when \x — y\ < T^/“, 
that is when the diagonal regime holds w.r.t. to the usual stable parabolic scaling. 
This is again the effect of the fractional in time derivative analyzed in small time, 
see (12.141) . Let us mention that there is some continuity w.r.t. to the stability index 
for the induced singularity, at least when d > 3 since the diffusive case provides a 
factor (T&\x — y\ d ~ 2 )~ l for that contribution in pp{T,x — y). 
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4. Proof of the Main Results. 

In the following we denote by c a generic positive constant that may change from 
line to line. Explicit dependencies for those constants c, mainly on assumptions 
(Ac) or (As), /3, are specified as well. 


4.1. Decomposition of the Error. For given points ( T,x,y) £ K.+ x (K. d ) 2 , we 
get from (12.141) and (13.21) that the error writes: 


£(T,x,y) 


(4.1) 


^ x z p 

[ Pzp (T, u)p(u, x , y)du - V] P [Z^' h = L]p h {U , x, y), 
JR+ 



Pzp ( T , u)(jp{u, x, y) - p(4>{u), x, y))du 


+{5Z P l Z T h = ti](p{ti,x,y) - P h (ti,x,y))}, 

i> 1 

£i(T,x,y) + £ 2 {T,x,y), 


denoting for all u £ (j>{u) := inf{t* := ih : ti > u}, where we have used 

(13.21) for the last but one equality. Let us specify that £\ is a term involving the 
time sensitivities of the density of the SDE (12. 1 2[) . which does not depend on the 
approximation scheme. The contribution £2 is associated with the spatial sensitivity 
involving the discretization error for the considered scheme. 

The idea is then to use some known asymptotics of the density S^’ + (see e.g. 
Zolotarev [Zol921 or Hahn et al. 11K1J111 ) which provide those for p Z f>. On the 
other hand, the quantities | p(ti,x, y)—p h {ti , x, y)\ have been thoroughly investigated 
in the literature for the diffusive case, see iKM0(4 . }GL081 . The results established 
therein give that for a general approximation scheme of type m the difference 
is controlled at order h 1 / 2 which is the convergence rate in the Gaussian LLT, see 
e.g. Petrov [ Pet05 _. If the scheme already involves Gaussian increments (usual 
Euler scheme of the diffusion) the control is improved and the convergence rate is 
h, see also jKM02l . In the strictly stable case, i.e. a < 2, the convergence rate 
of the Euler scheme, for which the exact stable increments are considered in m, 
has been investigated in |KM10] . The proof of Theorems 13.11 and 13.21 immediately 
follows from equation SID and the following Lemmas. 


Lemma 4.1 (Time Sensitivity in the Error). Under assumption (A/)) (Diffusive 
case) we have that there exists c := c(/3,(Ad)) > 1, s.t. for all T > 0 ,(x,y) £ 
(M d ) 2 , x^y: 


\£i(T,x,y)\ < ch 
exp (cT , ^/( 1+/3 )) 


exp(cT /3 / 2 ) ( c 1 \x — y\ 2 


TP\x — y\ d 


exp(—- 


TP 


rpPil+d/2) 


■ exp 1 — c 


„-i 


\x-y \ 2 

tp 


2^1 V(2-/3) 


f ( % 2 

\{TP/2\x-y\ 


Id>3 
\x - y\ 2 


pp(T,x-y) + j^pp{T,x 



< ch 
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Under Assumption (As) (Strictly stable case), there exists c := c(/3, (As)) s.t. 
for all T > 0, (x, y) £ (M d ) 2 , x ^ y: 


„)i < c 4 expfcr- 5 ") (fasqj +)+ 444 r /11 ,"“ (1 ~” ) 

1 v v ; \T^|x-y| d |x-y| d + Q J T^( 1+d / a ) V |x - y\ d+a 


< ch 

where w = — A 1. 


I \x-y\<TP/ a \x-y\>T?/ c 


+ t 7 y ) + ^w(r,x - y) 


Remark 4.2. Observe that in the stable case there is no loss of concentration as in 
the diffusive one. The second term of the stable bound in Lemma f4.11 corresponds 
to the usual estimate for the time derivative of the stable density, see Lemma [4.101 
below and [KolOO] . The first one corresponds to the additional spatial singularity 
induced by the fractional derivative in time in the diagonal regime \x — y\ < 
(stable parabolic scaling at time T^). 


We now give two Lemmas concerning the spatial sensitivity that involve the 
analysis of the discretization error and the time randomization. 


Lemma 4.3 (Spatial Sensitivity: Diffusive Case). Under assumption (A d) there 
exists c := c(/3, (A d)) s.t. for all T > 0, (x, y) £ (R d ) 2 , x y: 

- For the Euler scheme: 

(4.2) 

|£ 2 Eu 1 (T, x, y )| < ch{ | Id= ^ * j^(T, x - y) + x - y) J 

•— c/l£* j 5 ? Space(^'? % y) i 


where we write £|" ul , with a slight difference w.r.t. the definition in dUD, to specif¬ 
ically make the distinction with general Markov Chains. 

- For a general Markov Chain, we have that for all e £ (0,1/5) there exists c e := 
c £ ((A Dim )) s.t. for all T > 0, (x, y) £ (R d ) 2 , x/i/: 


\£- 2 (T,x,y)\ < c E 


h 1 ' 2 


(4.3) 


Id=l + Id>3 . Id=2 . 

\ X -y\ + J 'p /2 j qm -( d - 1 )+ Id =i’p( T ' X HI 

expfcrW 2 ) 


+ TPj2Vrn,p(T,X-y) ^ + 


exp l —c 




+ 1 


TP\x - y\ d ~ 2+Id ^ 

• — y\ \ — 2(m— l) + (d—2)+I <J < 2 


(/l 1 /^— e )l/ 2 


:— C e {/l ^ Bp, Space.LL t(^) X V) + ^"/3,Space,N oLLt(^ X V^)}- 


Lemma 4.4 (Spatial Sensitivity: Strictly Stable Case). Under assumption (A 5 ), 
there exists c := c(/3, (Ag))> 1 s.t. for all T > 0, (x,y) £ (R d ) 2 , x ^ y: 


|f 2 Eul (T,x,y)| < c/i( exp(cT^) 




TS \x — y\ d 


\x-y\>TP/° 

\x—y\d+ a 


exp (cT^/D-^O-g))) 
'J'Ptl + d/ a) \/|ai — y|d+a: 


< 


C ^| ^ 1 ^ | x—y\ a 7 - 7 X ~ V) + WPl 3 ^, X — y) j. 


Remark 4.5. We emphasize that in the strictly stable case, both time and spatial 
sensitivities yield the same error. This comes from the fact that the short time 
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singularity for the error expansion of the Euler scheme, exactly matches the time 
singularity of the derivative of the stable heat kernel, see equations (14.321) and 

(EH]). 

4.2. Asymptotics of the PDF of the Inverse Subordinator. As a preliminary 
remark, let us mention that in the particular case /3 = 2~ n , n € N*, for all T > 0, 
the random variable Z^, inverse at the level T of the subordinator S@’ + can be 
easily and explicitly simulated from the following Proposition. In particular, the 
auxiliary discrete approximation Z^' h is not needed in those cases, and the time 
sensitivity error (T : x,y) in (14.11) would vanish (considering a straightforward 
extension of the discretization scheme for the last time step). 

Proposition 4.6 (Identity in law for (3 = 2~ n , n € N*, T > 0). 

(4.4) Vn e N*, P = 2~ n , Z^ (1 = } \Z2\B 1 (V2\B 2 {V2\B 3 (• • • (V2\B n (T)\) ■ ■ ■ )|)|)|, 
where the (B l ) i& ^i... are independent Brownian motions. 

Proof. Let us start with n = 1. In this case, it is well known that, recalling that 

(law) 


t> 0 


(sr + ) t >o is the standard subordinator with index one has (^Sf 

( t_ t_ ) , where t_l_ := inf{s > 0 : B s = -4=}, stands for the hitting time of level 

V V2 /1 >0 C2 v2 

■^= for the standard Brownian motion B. Indeed, considering for a given A > 0 the 

exponential martingale (exp(-\/2A B s — As)) , the stopping theorem gives that, 

V / s>0 

for all s > 0, 

E[exp(-\/2Ai?sAT t ~ A (s A r_t_))] = 1 
Write now for allt > 0: 


==> E[exp(—Ar_f_)] = exp(— tX 1 ^ 2 ). 

S—>-+00 V2 


P[z| <t\ = P[Sf’ + >T} = P[r^ > T] = P[ sup B s < ~^=\ = F[V2\B T \ < t\, 

sG[0,T] V2 


y/2 


(law) 


using the well known Levy identity sup se r 0 T i B s = \Bt\, T > 0, for the last 

equality, see IRY99] . This proves (14.41) for n = 1. Introduce now a family of i.i.d. 

random processes ((>5i(t))t>o) i>1 with law (Sf' )t>o- It is then easily seen that 
S n (t ) := S n (S n -i(...Si(t))) is a standard stable subordinator of index /3 = -L. An 
immediate induction indeed gives: 

VA > 0, E[exp(-A5 n (S n _ 1 (...5 1 (t))))] = E[exp(-A 1 / 2 5„_ 1 (...5 1 (t)))] 

= exp(—A 1 ^ 2 h). 

Assume that identity (14.41) holds for a given neN* and let B n+1 be an additional 
independent Brownian motion independent of the (i?*)j e [i >n ]. Write: 


ns n +i(t) >T} = P [S n+1 (S n (t)) >T} = : 


sup < Ml] 

s6[0,T] v 2 


= F[C n (V2 sup B n+1 (s))<t], 

se[o ,T] 


(4.5) 
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where C n stands for the inverse of S n . Now, the induction hypothesis exactly gives 
that for all u > 0 

C n {u) °= } V2\B\V2\B 2 (V2\B\- ■ ■ (V2\B n (u)\) • • • )|)|)|. 

The proof then follows from the Levy identity and (14.51) since 5„+i is a subordinator 
of index l/2 n+1 . 

□ 


Remark 4.7. It can be shown similarly that for given (/S*)j£[ 1)n j £ (0,1)", n £ N*, 

if ((S'i(t)) t >o)are independent subordinators s.t for i £ [1, n], (S/(t))t>o ^ 

(S'f i ’ + ), then S n (t) := S n (S n _i(...Si(t )) is a subordinator of index /3 := ]/[" = i ft 
and 

Z T 

The delicate part is that this relation cannot be easily exploited from the practical 
viewpoint, except if /3j = i, Vi £ [1, nj. 

In particular, it follows from Proposition 14.61 that, for a given spatial Markov 
process {X t )t >o with generator L of the form (12.111) . the solution of equation (11.11) 
with [3 = 2~ n writes for a continuous function / and for all (T, x) £ R + x M. d : 


U(T,X) : E[f(X^-n)] ■ E [/("\^| B l(V2|B2(^| B 3(...(^| B n( T )|)...)|)|)|)]> 

generalizing Theorem 2.2 in Beghin and Orsingher p3O09l . 

We now state the Lemma concerning the required asymptotics for the laws of 

Zj rj~\ ^ Zj rj~t 


Lemma 4.8 (Bounds for the density of the inverse subordinator and its discrete 
approximation). For all T > 0, 

(4.6) Pzd (T, u) < 9(T,u), Vi £ N*,P [Z$ h = U\ < Cph9(T,ti), 


where: 



for some constants cp, Cp > 1. 


Remark 4.9. The above control is the generalization of the well known Gaussian 
bound for /3 = 1/2 recalled in Proposition [476] 


Proof. Let us first recall the following asymptotics for p s p,+ (l,-) from equation 
(2.2), (2.3) in IHKU11I . see also fZol92l : 

Ps«.+ (l, v) ~ exp (—[1 — /3 ](v//3)' S//(/3 “ 1 - ) ^ , v —> 0 + , 

(4.8) K = 1/\/ 27r/3(l — f3), 

p S 0 ,+ ( 1 , v) ~ v +°°- 












WEAK ERROR FOR FRACTIONAL IN TIME P(I)DES 


15 


Write now from (12.51) recalling the identity S'f ,+ > 0: 


P^{T,u) = -d»[f -±^ Psf> , + (l,J^ 


1 


(4.9) 


Pvt+W 

T 

Pvt+Vf* 


PSP’+ 1> 


Ps> 3 ’+ 
T 


(*’ J/p) dy+ 1 d v' PsB ' + (*’ Jjp) ydy 


V/3 


U 


integrating by parts and using (14.81) . 
Hence, 


P Z p (T, u ) 

Pzf< (T, u) 
u 


T~P 

r(i-/3) 


, u —> 0 H 


l 

T J 
+oo 


\ PI ( 2 ( 1 —/ 3 )) 


Ku^Tf exp ^-[1 - p\ (pf l{1 0) /( ^ . 


The above controls and elementary computations give the first bound in (14.61) . 
Indeed, there exists cq := cq(P) > 1 s.t., 


P z e(T,u) < 

P z e( T ’ u ) < 


T -0 / T 
c ni-py u - U 

^ s /3/(2(l-/3)) 


cK 


/ Q\P/W~P)) / f u xl/(l-/3)\ 



u > (c 0 Tf. 


For the last bound we observe that for P < L and u > (cqT)P, [u t—l ]/[T 2 ( 1 -' 3 ) ] < 
STi-fl= cT~P, with a constant c := c(P). For p > i, we obtain similarly 

that [u 1 ~p ]/[T 2 ^-^ ] < i -/ 3 . The last contribution can be absorbed by 

the exponential so that we are back to the previous case. For u € [(eg 1 T) /3 , (coT)^], 
the statement follows from (14.91) and the positivity of P s p+ (1, ■) on compact sets. 

From the control pze(T,u) < 6(T,u) and the relation (13.21) . we indeed derive 
that for all i £ N*, P[Z^’ h = ti} < Cph9(T,ti) for some Cp > 1. □ 


For the rest of the analysis we will need some controls on the p.d.f. of X t and 
their time derivatives in both the diffusive and strictly stable case. Such results 
can be found in Friedman |Fri64j under (Ad) or [KolOOJ under (As). For the 
reader’s convenience we also recall in Appendix[A]some points about the parametrix 
expansion from which these controls are derived. Let us emphasize that, for our 
analysis, we need to specificy explicitly the dependence of the constants in time 
since we are led to integrate over arbitrarily large time intervals in (12.141) . 

Lemma 4.10 (Non degenerate heat kernel bounds). Under Assumption (Ad) (dif¬ 
fusive case), there exists (c,C) := (c,C)((Ad),cI) > 1 s.t. for all (u,x,y) € 
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2. 


p(u,x,y) < Cexp(Cu 1/2 )g c (u,x - y), 

ip , m ^ exp (Cu 1 / 2 ) 

|d„p(u,x, 2 /)| < C- g c [u,x-y ), 


(4.10) 


g c (u,x-y) := 


1 , lx — y | 2 

• exp(— o ). 


(2-KCu) d / 2 2cu 

Under Assumption (As) (strictly stable case), there exists C := C((As),d) > 
1, w := ^ A 1 s.t. for all ( u,x,y ) £ x (R d ) 2 : 

p(u,x,y) < C exp(Cu u )ps(u,x - y), 


d u p{u,x,y) < c eXP ^ Cu ^ p S {u,x - y), 


(4.11) 


p s (u,x-y) = 


,.d/a 


i + 


\x - y | 

v l/cx 


— ( d+c 


where the constant c Q is chosen to have f Rd ps(u, z)dz = 1 . 

Observe that the above controls are exponentially explosive in time. 

We now give a fundamental lemma for the error analysis. The quantities below 
derive from the time sensitivities or the previous error bounds from [KMOOj , |KM02j 
used to control the terms £\, £i in (14.11) respectively. 

Lemma 4.11 (Fundamental Bounds). Let T > 0 be given. The following controls 
hold: 

- Diffusive case: There exists c := c((A/j)) s.t. for 7 £ {0,^,1} and for all 

(x,y) £ (R d ) 2 , x^y: 

r+00 

/ u~ 1 pzP ( T , u)g c (u, x — y) exp (Cu 1//2 )du 

Jo 

+ X! *r 7p [^T ' h = ti]9c(ti, x-y) exp (Ct] /2 ) 


< c 


{<■ 


i> 1 


Id=2 


+ 7- , 1 ^ 1 x - y) 

/r* — Qi\u-'Y>O r r 9. J W=1 / 


Irf>3 _!_ 

2 - — y? 1 \x — — y\h>oTTi I r 


(4.12) 




For the Markov Chain approximation, for all e £ (0, |), there exists c e := 
c E ((A D ,m)) s.t. denoting for all (x,y) £ (R d ) 2 , 

\ —m 

-y\\ 


q m (U,x-y) : =-dH \ 1 + — 


1/2 


where c m is s.t. f Rd q m (ti, z)dz = 1 , we have for all x ^ y: 

T. q m (ti,x- y)t~ 1/2 exp(Ct] /2 )P[Z^ h = f;] 


< c E 


i> [ft. -4 / 5-6 ] 

Id=l + Id>3 , Id=2 ) 


1 


_ y| + rpp /2 j < l™-[(d-i)+i d =i],p(T,x y) + Tpi 2 q mt p(T,x y) 

:= c e^/3fs P ace,LLT(^ x ~ V)■ 


(4.13) 
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- Strictly Stable case: There exists c := c((As)) s.t. for 7 £ {0,1} and for all 
(x,y) £ (R d ) 2 , x^y: 


(4.14) 


r+oo 


/ u 7 p Z / 3 (T,w)ps(u,z - y)exp(Cu UJ )du 
Jo 

+ ^2tpF[Z^ h = ti)p s (ti,x- y ) exp(C^) 


i> 1 


< C 


{( 


[ \x-y\<T?/°‘ ^\x—y\>TP/ a \ „ , s . 1 - (rr A 

+-77^-Jp i 8 (r } a:-i/) + ^p j 8 (r,a:-y)|. 


As a direct corollary of equation (12.141) and Lemmas 14.101 and 14.111 taking 7 = 0 
we obtain some controls on the density of X z p. Precisely, we have: 


Corollary 4.1 (Some density estimates for the spatial motion along the inverse 
subordinator). There exists c := c((A),/3, a) > 1 s.t. for any fixed T > 0 the 
following estimates hold: 


Px a (x, y) < c(p 0 + pp)(T , x-y), 

Z T 

in both the diffusive and strictly stable case with the associated definitions of pp,pp 
in Theorems 13.11 and 13.21 


Remark 4.12. In the above bounds, from the definitions oip^,p^ 1 we have in some 
sense the usual scaling in time for the Gaussian and stable heat-kernels at time 
T' 3 . We indeed observe that the dichotomy between the diagonal and off-diagonal 
regime is performed depending on \x — y\ < K(TP) or \x — y\ > A'(T /3 ) 1 /“ 
for some positive constant K. In the diagonal regime we also have, as indicated 
previously, an additional spatial singularity whereas in the off-diagonal one, the 
concentration is modified in the diffusive case due to the integration of the density 
Pzt>(T,u ) for large u in (12.141) . 

In the Gaussian case, those bounds are coherent with those of [EK04] and our 
approach can be viewed as an alternative probabilistic strategy to derive estimates 
for the fundamental solution of fractional in time diffusive PDEs. For the strictly 
stable case, to the best of our knowledge these results seem to be new. Some lower 
bounds are derived in Appendix [Bj 

Proof of Lemma \f.ll\ 


- Diffusive case. We focus on the first integral in the l.h.s. of (14.121) since the 
other term can be handled similarly (Riemann sum approximation). Exploiting the 
bounds from Lemma 14.81 we split the time integration domain at the characteristic 
time scale and write: 


r+oo 

(4.15) Q Jt p(T,x-y):= u^pze (T, u)g c {u, x - y) exp(Cu 1/2 )du 

Jo 

csexp(GT^/ 2 ) f c~ 1 \x — y\ 2 \ f T P du ( c -1 la: — y\ 2 \ 


+ 


2L 

tp 



exp(Cu 1,/2 ) exp 



:= < 2 7 ,/ 3 , s (T, x-y) + Q lt /3,i{T, x - y), 
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where the subscripts s, l respectively denote the contributions in short and long 
time. For the contribution Q lt p tS (T, x—y), we have to analyze the integral I (T^ 3 , 7 , \x— 

y\,d ) := / 0 T u ^+d /2 ex P( —£ 4 ~ )• Observe that except if d = 1,7 = 0 there is 

a singularity in small time in the integral which needs to be compensated by the 
spatial component \x — y\ ^ 0 in the exponential. This yields the spatial (diagonal) 
singularity. 

To equilibrate the time singularity we will thoroughly rely on the fact that the 
density of the hitting time of level a > 0 for the scalar Brownian motion writes (see 

e.g. fTTvTTm i: 


(4.16) 


Vit > 0, fa{u) := 


2 a 


( 27 TO 3 ) 1 / 2 

If 7 + d/2 > 3/2 we write for all \x — y\ ^ 0: 


6XP 1 ~2u 


i(T 0 nA x -y\,d) 


j'& 


< 


Jo 

c 


du\x — y\ \x — y\ d+21 3 

u (d+ 2 7 -3)/2 


exp 


y I s 


\x — y\ d + 2 (j- 1 ) 


’0 


du\x — y | 

(2 tto 3 ) 1 /2 


exp —c 


jr-l 


4 u 
r-y\ 2 


1 


2 it 


< 


|x — y | d + 2 7 —2 
c 

|x — y|d+2(7-l) ' 


■ If 7 + d/2 < 3/2 which actually happens for d = 1,7 £ {0, or d = 2 and 7 = 0 
we write: 

J(T^, 0 , |x — 2 /|, 1 ) < cT^/ 2 , 


/(T^, 0 ,|* - 2 /|, 2 ) = /(T^, 1 / 2 , |a; — y\, 1 ) < 


rp/3/2 


2/1 


du|a; — y | 
7/3/2 


■ exp 


k - y | 5 


< 


cT ^/2 

I*- 2 / 1 ' 


The previous results can be summarized as follows. There exists c := c((Ad),/3) 


s.t. for all (T, x, y) £ x 


- 2/1 7^0: 


(4.17) 


Qj,p,s(T,x- y) < cexp(cT /3/2 )exp 1 ^ 2T f ^ 

}• 


"Id=l,7=0 I<J=2,7=0 + Id=l,7=1/2 


Id+27>3 


l T/3/2 T^/ 2 |a; — 7/| T/3|a; - 7/| d +20-i) 

The restriction \x — y\ 7/ 0 is not necessary if d = 1, 7 = 0 . 

On the other hand, from Young’s inequality, there exists cp > 0 s.t. for all 
u, e > 0 : 


u 1 ' 2 = 


(e/i ) 1,2 | + e -./(i+/» T />/(. + « J . 

We thus get from (14.151) that, for e sufficiently small, there exists cp s.t.: 

Q^i(T,x- y ) < T0[1 C ^ +d/ 2) exp(^T^( 1+ ^)) 


/■+00 

x / du exp ( —c 

Jtp 


-1 k~ 2 /P 
2 u 


exp -- 


Ca 1 r U '|1/(1-/5) N 


I—V 

2 It/ 3 / 
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Now two cases occur according to the dichotomy indicated in Remark 14.121 If \x — 
y \ < T 0 / 2 then setting u := u/TP we get Q 7 > ^ i i(T, x-y) < TP ^+ d / 2) exp(c / 3 T /3 /( 1+/3 )). 
If now \x — y\ > T^/ 2 , we see that the two exponential contributions in the integral 
are equivalent if Ki{\x — i/\ 2 tP/ ( - 1 ~P^} 2= ^ <u< K 2 {\x — y\ 2 TP^ 1 ~P^} 2= ^, 1 < K\ < 
K 2 < +oo, for which 



On the other hand, for u G [T&, Ki{\x — y\ 2 TP^ x ^} 2 -?] we have: 


exp l — c 


„-i \x-yf 


2 u 


< exp ( — ( 2 c«i) 


■~y \ 2 

tp 


1/(2-/3) 


For u £ [m 2 {|* — y\ 2 TP/^ P^} 2 -?,+oo), it suffices to bound 


exp 


c « ( u 


4 {^} 


V(l-/3) 


< exp — 


c /3 1 i/(i-/3) 


2 \ l/(2-/3) 


T/ 3 


Thus, setting as in the diagonal case u = u/TP, we derive that there exists 
c i g((Afl),Ki,K 2 ) > 1 , s.t. for all ( T,x,y ) £ R+ x (R d ) 2 : 

(4.18) 


(T, •£ 


2/) ^ T /3(7^/2) exp(c^r^( 1+ ^)exp ^-Cg 1 


>- v\ 2 ' 
tp 


1 /( 2 -/?) 


C/3 := 


The statement (14.121) follows from (14.181) and (14.171) . 

- Markov Chain. 

Setting ic := |~/i _4 / 5_e ], j := [T^/h], let us write from Lemma FOsl 


Qi,p(T’X~y) '■= T p i Z T h =ii]exp (Ct 1 / 2 )^ 1/2 q m {ti,x — y) 

i>ic 


< 


J 7 ?(l+(d+l)/ 2 ) 

(4.19) 


exp (cTP/V+P')) J 


TP 

+oo 


exp {OTP/ 2 ) / “ u~( d+1)/2 ( 1 + 


\x-y\ 


1 + 


- y\ 

7,1/2 


exp( 




u 1 / 2 

l/(l-/3) 


drt 

)du 


: = Q±,p,s( T ’ x - y) + Qi,p,i( T > x ~ y)> 


using as above Young’s inequality to get the upper bound and replacing, up to 
modifications of the constants, the Riemann sums by integrals. 
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We now derive setting t = u/\x — y \ 2 and up to a modification of c: 




< 


TV 3 


exp (cT^ 2 > 


t. 


T^d=1 


X — y|d-l+Id=i 
bld=l 


ti/\x-y\ 2 / i 

(d+l+Id = i)/2 ( X + 

ic/N-yl 2 


fl /2 


di 


t? 


— 7^0 ex P(cT^^ 2 )-- ~ - 

TP \x — y\ d 1 + 1 < i = 1 


/•ti_/\x—y\ 2 / I \ — m+(d— 2+I<; = i) 


dt 


T P/\ 


■d=l 


' - 2/1 


'ti c /|x- 3 /| 2 

Id>2 
' - y \ d ~ 1 


X 1 + 


x-y\ 
TPt 2 


-m+(d- 1)+Id=i 


< ^ exp (cTP/ 2 ) 

(4.20) 

On the other hand if \x—y\ < T ^/ 2 (diagonal regime) we readily get Q i p t i(T, x— 
y) < ce ^ 0 (d+i )/2 ^ ■ If now \x — y\ > T^/ 2 , splitting as above for u G [f,, Ki{|x — 
y^T^/C 1- ^)} 2 -^] and u > Ki{|a; — ypT^ 1- ^} 2 ^ yields : 

n l/(2-/3)\ ~ m 


QbdA T ’ X ~V)^ T p(d+ 1)/2 exp(cT^ 1+ «) ^1 + 


k ~ 2/1 

T p/2 


(4.21) 

Equation (14.131) follows putting (14.211) and (14.201) in (14.191) . 

- Strictly Stable Case. We again focus on the first term in the l.h.s. of (14.141) . 
The discrete sum can be handled similarly. Analogously to (14.151) . replacing g c by 
ps and the exponent 1/2 in the exponential by u>, we write from Lemmas 14.81 and 

mm 

cpc a exp (CTP U ) f T 13 du 1 


Qf,p, s (T,x - y) < 

Two cases are to be distinguished. 


TP 


]o u^+ d / a (1 + \ x ~y\ )d+ 


-Ifd = l ,7 = 0 ,a>l then the above diagonal singularity is integrable and one 
gets: 

i pT 13 


Qo,/3,s (-1/ 2/) A 


cpc a exp (CTP U ) 


du 


TP 


(4.22) 


< 


(1 + ^) 1+a do « 1/0 

cc a exp(cT /3w ) 1 


TNa (l + ^) 1+a ' 

In particular, this estimate holds even for x = y. There is no spatial singularity. 

- For all the other cases, i.e. d = 1 ,7 = 1, a > 1 or d = 1, a < 1, d > 2 for 7 £ {0,1}, 
we need to consider i/y in order to equilibrate the time singularity. In small time 
we still distinguish the diagonal and off-diagonal regimes. Precisely, 

- If |cc — z/| < T^/ a , then: 
cpc a exp (CTP U ) 


( 1 r \ x - v \°- r 

{W^W^L duu + L 


Qj, 0 ,s(T, x — y) < 


TP 


\ x -y\ d +a 

(4,i3) - C/3Cexp ( cT/3 ^{tP\x - + TPh+ d / a ) }■ 


rp/3 

\x-y\ c 


du 
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- If \x - y\ > TP/° 


x y ) ^ 


cpc a exp(CT^") 


1 




tp 


(4.24) 


< c / 3 cexp(cT /3w ) 


\x - y\ d+a 
j'Ptp-'i) 


duu 1 7 


- y I 


d-\-ct ’ 


Let us mention that, in the diagonal regime, we get as in the diffusive case an 
additional spatial singularity. 

Let us now deal with Q lt p t i(T,x — y). We first write from Young’s inequality 
that there exists cp s.t. for all u,e > 0: 

u “ = ( £ ^r {£ ~ lTP)U - ^ {( e ^) 1/(1_/3) + (e-Wr/V-* 1 -™ |. 


Thus, taking e small enough yields: 


x - y) < c 


exp(cT /3aj /( 1 - aj ( 1 -^ 3 ))) [+° 


x exp 


TP 

2 XtP ) 


Jtp 

i/(i-/?)' 


dzt(-j 


- 2/1 




A z t - d/3 /“)u- 7 


(4.25) 


< cexp(cT /3 “/( 1 -" J ( 1 - /3 ») f --A T^A T _/37 . 

VF- 2/1 +a / 


|x - z/| d +“ 

The result follows from (14.221) . (14.231) . (14.241) . and (14.251) . 


□ 


4.3. Time Sensitivity: Proof of Lemma 14.11 The lemma, which gives the 
control on the time sensitivity term £\(T,x,y) in the decomposition error (14.11) . 
is proved using the bounds of Lemma 14.111 We simply write from Lemma 14.101 
and Lemma 14.81 for the first time step, that in the diffusive case, there exists c := 
c((A d )) s.t. : 


< c 


hexpiCh 1 / 2 ) 


exp l —c 


|£i {T,x,y)\ 

-i \ x ~ y \ 2 


TP\x — y\ d 

r+oo 

+h / Pzp (T, u)u~ l exp(Cu 1 ^ 2 )g c (u, x — y)du 
Jh 


whereas in the strictly stable case, there exists c := c((Ag)) s.t. : 


|£i( T , x,y)\ < c| + h j h Pzf>(T,u)u 1 exp(Cu u )p s (u,x - y)dv)j . 

Recalling that in the strictly stable case we assumed h 1 ^ < T and h l / a < \x — y\, 
the result then follows from the above controls and equations (14.121) . (14.141) in 
Lemma 14.111 


4.4. Sensitivity in Space: Proof of Lemmas 14.31 and 14.41 
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4.4.1. The case of a diffusive spatial motion. The key control for the analysis of 
the term £ 2 (T,x,y) in ( 14.111 is provided by the following Lemma. 

Lemma 4.13 (Control for the Spatial Error). For our current approximation schemes 
in HO we have: 

- Under (A/yEui), i- e - f° r the Euler scheme, there exists (C,c) := (C,c)((Au),d) > 
1 s.t. for a given h > 0 and for all i £ N*: 

\(p-PEu\)(.ti,x,y)\ < Ch {~Tp. V exp(Ct] /2 )}g c (U,x - y), 

(4.26) tl 2 

s “ (i ^ ) = (2dy^ exp ('S")' 

standing for the usual d-dimensional Gaussian density with variance cU- 

- Under (Ac im ), i.e. general Markov Chain approximation, there exists C := 
C((A Dtm ),d) > 1 s.t. if ti > h s , S < 1/5: 


(4.27) 


\(p~P h )(ti,x,y)\ < Ch 1/2 {^ V exp(Ct] /2 )}q m (ti,x 


0 < q m (ti,x- y) := c m t. 


-d/2 



e J 


y). 


Let us mention that q m enjoys the same parabolic scaling as g c but has poly¬ 
nomial decay. We observe the usual convergence rates in h and h 1/l2 respectively 
for the Euler scheme and the Markov Chain approximation. For general Markov 
Chains, this is the rate of the Gaussian LLT (see [Pet05l . [ KM0Q i. [KM09j h For the 
Euler scheme, the innovations are already Gaussian so that the first term disappears 
in the previously mentioned LLT, yielding a contribution of order h. 

The important point of the previous Lemma is that it specifies the behaviour 
of the constants in short and long time. This allows us to balance, as for the 
time sensitivity, those constants with the controls on the laws of Zj,, Z^ h given in 
Lemma ITT51 thanks to Lemma T4. Ill 


Proof. For the Euler scheme, the short time behaviour of the error has already been 
established in C1.08 (see Theorem 2.3 therein) whereas the exponential bounds 
can once again be derived from the explicit form of the parametrix series used 
to investigate the error in [ KM021 . For the sake of completeness, we recall these 
aspects in Appendix lA.il 

For a general Markov Chain approaching the diffusion, i.e. for which the inno¬ 
vations are not necessarily Gaussian, to enter a Gaussian asymptotics specified by 
the Gaussian LLT, a certain number of time steps is needed. This is why we impose 
the condition U > h s , S < 1/5 which is the one required in IKM09] (see assumption 
(B2) therein) to establish the Gaussian LLT in short time. The behaviour of the 
constants in large time can also be derived from jKMOO] (similarly to the procedure 
presented in Appendix IA.1I) . For the control in short time we refer to Remark 1 
after Theorem 1 in [KM091 . Anyhow, for the reader’s convenience, we recall in 
Appendix IA.1I the approach developed therein. □ 
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We thus directly get from (14.111 and (14.2611 in Lemma 14.131 that there exists 
(c, C) := c((Ad))> 1 s.t. for all T > 0, (x, y) £ (K d ) 2 , x ^ y: 

(4.28) \ £ l P ul ( T > x -> y)l < Ch E p i Z T h = **] exp(Ct■' /2 )t~ 1/2 g c {ti, x-y). 

i> 1 

The result (14.21) for the Euler scheme now follows from equation (14.121) in Lemma 

sm 

For the Markov Chain approximation, setting ic '■= [A _4 / 5_e ~| we write using 

(TPD) : 


(4.29) 


\£ 2 {T,x,y)\ < Ci Y W[ Z T = ti)\(p-p h )(U,x,y)\ 

+ VhY p [ Z T h = exp {Ct\ /2 )-^q m (ti,x - y) j 

4 4 t/i ' 


i>ic 


< c e {£ 2 i(T, x, y) + ^£^ SpaC e, LLT (r, x - y)}, 


exploiting equation (14.131) in Lemma 14.111 for the last inequality in the previous 
r.h.s. It therefore remains to control £21 (T, x, y). 

For this contribution, we do not have comparison results between the two den¬ 
sities. The technique thus consists in controlling each of the terms. Observe to 
this end that a simple parametrix expansion for the scheme, similar to the one per¬ 
formed in (KMOOl using Lemmas 3.6 and 3.11 therein, (see also Appendix lAll. would 
give that there exists c := c((A/y m ), d) > 1 s.t. for all i £ [l,«c]> (x,y) £ ^ d ^ 2 - 


(4.30, (l + t#) 


— 2 (m— 1 ) 


We thus derive from the definition in (14.291) and equations (14.61) . (14.101) . (14.301) that: 


£21 {T,x,y) < ft exp(cT^/ 2 ) 


E* 

i=1 


-d/2 


exp - 


\ x ~y? 


ti 


+ 1 


\x-y\ 


— 2 (m— 1 ) 


t- 


1/2 


From the arguments in the proof of Lemma 14. Ill we can write for all d > 1: 

,|ld<2 


t? 


f 2i (T,x,y) < ^exp(cT^ 2 ) |a; _^_ 2+I ^ 2 


exp 1 —c 


\x-y\‘ 


t 


+ 1 


(4.31) 


<£ 


M.e 


— /3,Space,NoLLT 


lc 

(T,x-y,h), 


\ x - y I 

t 1/2 

L ic 


-2(m-l)+(d-2)+I d < 2 ' 


recalling t lc = h 1 / 5 e . Putting together estimates (14.311) into (14.291) for the Markov 
Chain approximation concludes the proof of Lemma 14.31 
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4.4.2. The case of strictly stable driven SDE. Similarly to the diffusive case, the 
crucial point is the following Lemma. 


Lemma 4.14. Under (Ag) there exists C := C((Ag)), s.t. for a given U := ih we 
have: 


\(P~PEu\)(.ti,x,y)\ < Ch{— V exp (Ctf)}p s {t i: x - y), 


(4.32) 


Ps { U , z ) := 


d/a 


i + 4^ 

t' 


d+a. 


, Vz e R d , f ps ( U , z)dz = 1. 
Jw d 


Proof. The exponential control is derived similarly to the diffusive case, and comes 
from the specific form of the parametrix expansion used to analyze the error. The 
control in small time can also be derived from this representation, similarly to what 
occurs in |KM09l . To make this last point clear we give some details in Appendix 

m □ 


Lemma 14.41 now follows from Lemma 14.141 and equation (14.141) in Lemma 14.111 


5. Perspectives 


We did not consider in this work the case of a general Markov Chain approxi¬ 
mating stable driven SDEs. To do so the first step would consist in establishing a 
LLT for sums of i.i.d random variables belonging to the domain of attraction of a 
stable law. Such results have been thoroughly studied in the scalar case, see e.g. 
Mitalauskas and Statuljavicjus [MS76I or the monograph by Christoph and Wolf 
|CW92| . We also refer to the PhD. of Squartini |Squ07| and to [MPS07] (Cauchy 
case), for recent developments. We believe those results can be extended to the 
multidimensional case and would also transmit to a Markov Chain approximation 
through a continuity technique like the parametrix. Now such limit theorems could 
also allow to consider more general fractional like time derivatives that would be 
associated with the inverse of inhomogeneous stable subordinators (S^~)t>o with 
generators 

L t f>{s) = f ((f(s + u)-(j){s))g(t,u)-^. 

JR^ U ^ 

Eventually, an analysis of generalized fractional like derivatives involving as well the 
spatial variable like in Kolokoltsov’s work, see e.g. Proposition 4.4 in IKol09j . would 
require to establish LLT for approximations of the couple (.S) + , X t ) with generator 

/* du 

£t4>(s,x) = / {(j>(s + u,x)-(j)(s,x))g(t,u,x)—rr- 

J r ; u 1+a 

f dz 

+ / x + z )- 4 >{s, x)) f(t, x , z ) . 

JR d l~l 

This case is much more delicate in the sense that the complete coupling breaks 
the independence in (12. 13[) between the spatial motion and the process associated 
with the fractional like time derivative. We anyhow believe that the arguments of 
[Kol09j can be adapted provided the LLT holds. This will concern further research. 
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Appendix A. Derivation of the Error Bounds in Small and Large 

Time. 

We briefly explain in this section the bounds appearing in (14.261) . (14.271) and 
(14.321) that are crucial for our analysis in order to balance the indicated explosive 
behaviours with the decays of the density of the inverse subordinator, see Lemma 
14.81 Actually, the only control that needs to be fully justified is (14.321) in small time. 
The other bounds in small time are already established in the previously quoted 
papers. Also, the exponential bounds are in some sense classical. 

Let us consider the case of the Euler scheme associated with (12.121) first. The 
crucial point is that the densities p(f,, x, -),Peui(^u •)> °f respectively X ti in (12.121) 
and its Euler scheme Xjf in in starting at x, enjoy a parametrix expansion, see 
again |Fri64| , [MS 6 7] , [KMOOj , |KM02j for the diffusive case and [KolOO] and |KM10] 
for the stable one. We proceed with the simplest example of a non degenerate 
diffusion, but the results for the other cases can be derived similarly. Consider the 
dynamics (12.121) under (A^). The p.d.f p(t,x ,.) of Xf exists for every t > 0 and 
writes: 

(A.l) p(t,x,y) = p(t, x, y) H^{t,x,y), 

i> 1 

where for all ( t , x,y ) € x R 2d , H(t,x,y) = (L — L)p(t,x,y) where p ( t , x , y ) 
stands for the p.d.f. at point y of the frozen Gaussian or stable process Xf’ v = x + 
b(y)t + cr(y)Yt, with generator L. Note that the frozen density is always considered 
at the freezing point. In (IA.1I) for two integrable functions /, g : x (K. d ) 2 —> R., 

we denote / <g) g(t,x,y) := fgdu f Rd f(u, x, z)g(t — u,z,y)dz and also for all % > 
1 , := HO” 1 ) f® HM = f. 

A.l. Diffusion Case. It is known (see }Fri64l . (KM02j ) that, under (Ad), there 
exist (ci, C 2 ) := (ci, C 2 )((Ajj))> 1 s.t. for all * > l,t > 0, 

(A.2) \p® H {l) (t,x,y)\ < r [ Clt \ /9 s g C3 {t,x- y). 

r (* — 1/2) 

From equations (IA.2I) and (IA.1I) . it is easily seen that there exists (ci,C 2 ) := 
(ci,c 2 )((A£>)) s.t. 

(A.3) pit , x, y) < ci exp(cit 1/2 )c/ C2 (f , x-y). 

This explains the exponential control. This approach is due to McKean and Singer 
1MS67 . It can be easily transposed to the approximation schemes. The p.d.f. of 
Euler approximations enjoy similar properties as the one of the initial SDE (see the 
quoted references). 

The control in small time comes from the explicit form of the error decomposition 
which is similar to (IA.1I) . Namely, one has (see [KM02| in the considered case) 

(A.4) (p-PEu\)(ti,x,y) = ^(p® {L 2 — L 2 )p) (U,x,y) + h 2 R{U, x, y), 

where the remainder term satisfies as well |i?(tj, x, y) \ < Ci exp(cit,[' /2 )t !; X ^ 2 g C2 {ti,x— 
y). Also, we denote for i G {1,2 },L\<j>{x) := {L\<j>{x ))| 5=x , L^(j>{x) = (6(0, V0(o:))+ 
|Tr {a(^)D^.<j>(x)). Observe that L<j>{x) = L*<t>{x), but more generally the operators 
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do not coincide anymore when iterated. Precisely considering d = 1 to alleviate the 
notations, we get for all t £ (0,fi): 

(L 2 - Ll)p{t,x,y) = 

{b(x)b'(x) + -a(x)b"(x)}d x p(t,x,y) 

+{^b{x)a'(x) + a{x)b'{x) + ^a(x)a" (x)}d 2 p(t,x,y) 

(A.5) + ^a{x)a'{x)d^p{t,x,y). 

Now, the controls on the density and its derivatives imply under the smoothness 
assumptions in (A^) similarly to (IA.3I) that there exists c := c((A^)) s.t. for all 
multi-indexes a, / 3 , |a| + |/ 3 | < 3 for all t > 0: 

(A.6) \d%dPdp(t,x,y)\ < |ct[ C +|g , exp {ct 1/2 )g C2 {t,x - y). 

t 2 

We refer to Friedman |Fri64l for details. Hence, the most singular term in (IA.5I) is 
the last one. Precisely, 

(A.7) \a{x)a'(x)d*p(t,x,y)\ < ^exp {ct 1/2 )g C2 (t,x - y). 

Now, recalling from (IA.4I) that this control needs to be plugged in the time-space 
convolution p (§) [{Li 2 - Ll)\p(ti, x, y) = f 0 ' ds f Rd p{s, x , z){L 2 - L\)p{U - s, 2 , y)dz , 
two cases can occur: if s £ [0, £;/2] the control in (IA.7I) is not singular and gives the 
stated bound of t,- ' once integrated in time. On the other hand when s £ [fj/2, tj, 
some integration by parts need to be performed for the singular term yielding 
a third order spatial derivative in 2 of p{s,x,z). Thanks to (IA.6I) . which gives 
|<9?p(s, x, z) I < exp(cs 1 / 2 )^, (s, x — z), and the integration in time, we get again 

a time singularity in t i 1/ ” 2 for the main contribution in (IA.4I) . The remainder can 
be handled similarly (see [KM021 1. This analysis can be performed as well for the 
Markov Chain approximation, see [KM09I . 


A.2. Strictly Stable Case. The expansion (IA.1I) also holds under (Ag), see 
(KolOO] from which we have 

\p® H {l \t,x,y)\ < ^. Cl ^ p s {t,x-y), u:= f-Al) £ (1/2,1]. 

r(*-w) ) 

This therefore gives the indicated exponential bound p{t, x , y) < cexp{ct ul )ps{t,x — 
y) for c := c((As)). Now, let us denote by / the spherical density of the measure 
v (see assumption (ND)). We can rewrite: 


Lip(x) = L*^(x) = {b{x), V x ip(x)} 

dz 


+ / {v{ x + a{x)z) — ip{x) — (V x ip(x), cr(x)z)I| (T ( x)z |< 1 }/ 

Js. d 

= {b(x), v x ip(x)) + / {<p{x + z) - tp(x) - {S7 x ^p{x),z)l\ z \< 1 }Q{x,z)dz, 


where we denoted for all £ £ 


(A.8) 


0 (C>*) : =t^t 


r / v-'lOz 


cr _1 (C)^| d+a det(cr(C)) " 
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With the notations of the previous section (considering a < 1 to avoid the cut-off 
and alleviate the notations) we get: 

L 2 p(t, x, y) = {b(x), V x b(x)V x p{t , x, y)) + (b(x), D 2 x p{t, x, y)b(x)) 

+ {b(x), / ( V x p(t , x + z, y) - \7 x p(t, x , y))Q(x, z)dz) 

Jg. d 

+{b(x), / ( P(t, x + z,y) — p(t, x, y))V x 0(x, z)dz) 

jR d 

/ \[ {p{t,x + z'+ z,y)-p(t,x + z',y)}Q(x + z',z)dz] 

JR d L J R d 

-[ / {p(t,x + z,y) - p(t,x,y)}Q(x, z)dz] l©(a :,z’)dz’, 

jR d J 

(b(x),D 2 p(t,x,y)b ({)) 


+ 


L 2 p(t, x, y) = 


+{b(x), / (V x p(t, x + z,y) — S7 x p(t, x, y))0(£, z)dz) 
jR d 

+ [ \[\ {p{t , x + z' + z,y)~ p(t, x + z', j/)}0(£, z)dz] 

jR d l J R d 

- [ [ {p(t,x + z,y) -p(t,x,y)}e(x,z)dz]\<d(t,z')dz' 
JR d ) 


Thus, the difference writes: 

(L 2 — L 2 )p(t,x,y) = (b(x),S7 x b{x)V x p(t,x,y)) 


+(b(x), / (p(t, x + z,y ) - p(t, x, y))\7 x Q(x, z)dz) 

jR d 

[ {[/ {p{t,x + z' + z,y)-p(t,x + z’,y)} 
jR d l J R d 


+ 

x{0(x + z 1 , z) — 0(x, z)}dz] x 0(x, z')dz'. 

We mention that, even though we assumed b = 0 for a < 1 we kept the explicit 
dependence on b in the above computations for the sake of completeness. Let us 
now focus on the last term in the above equation: 

D(t, x, y) := / / {p(t,x + z + z',y)-p(t,x + z',y) 

JR d l JR d 

x{0(x + z 1 , z) — 0(x, z)}dz | x 0(:r, z')dz ', 

which is the only contribution if a < 1 and which would be the most singular for 
a > 1. The smoothness and non-degeneracy assumptions (S), (UE), (ND) on the 
coefficients yield: 

l V * e (*’*)l ^ |Tp+^ 

\D(t,x,y)\ < cl I \p(t,x + z + z',y)-p(t,x + z',y)\ 

J R d JR d 

I P'|>it 1 /“ + \ z '\^\z'\<\tV a 1 


I d+a 


\ z '\ d +a 


dz'dz , 
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where we have truncated with respect to the characteristic time-scale t 1 /“ (up to a 
constant) in the variable z'. Indeed, if \z'\ < \t 1 ! a we perform a Taylor expansion 
and exploit the previous bound on V a -0 whereas if \z'\ > \t x ! a we simply use the 
uniform bound of 0 deriving from its definition in (I A. 81) and the non degeneracy 
assumptions. Write now: 


\D(t,x,y)\ < c / \p(t,x +z +z',y) - p(t,x + z 1 ,y)\ 

•Ilz'lcit 1 /" JR d 


l\z'\< 

X , ! Z '\ . X 7— ■rrr~d Z r dz 


+c / \p(t,x + z + z',y) -p(t,x + z',y)\ 

J\z'\>kt x / a Jm d 

1 


-dz’dz \= {Di + D 2 )(t,x,y). 


Let us treat those two terms separately. Recalling from [KolOOj and the above 
bounds that \\7 x p(t,x,y)\ < ce *i/)f ^ ps(t,x — y) we derive: 

Di(t,x,y) < ^ [ p s {t,x + z' + 9z-y)\z\ 

x ^ I x 1 d J dz 

U/ld+a | 7 |d+ci az 




I p(t, X + z + z',y) - pit , x + z', y)\ 


\ z ' I 1 

x TZm^ dz ' dz ~ {Dn+Di 2 ){t,x,y), 


for some 9 := 9(x, z, z ', y) £ [0,1] in D n . In that contribution, since both z, z' are 
small w.r.t. the characteristic time f 1 /", we have: 


Ps{t,x + z' + Oz - y) < 


1 


< 


1 


(A.9) 

Hence: 

Du(t,x,y) < 


fd/a ^ \x+z'+6z-y\ y :+n ~ f-d/a (1 |x-y| y +n 

< c2 d+a p s (t,x-y). 


c exp{ct u 
t 1 / 01 


- ps(t,x - y) ( f - 


dr \ cexp(cf^) . . 

1 < -v-b— ~Ps(t, x — y). 


a / — XI-X 


For the contribution | p{t, x + z + z',y) — p(t, x + z ', y)\ in Di 2 {t, x, y), since \z\ can 
be large, we write directly, recalling as well that \z'\ is small and proceeding as in 

(ESI): 


| p{t,x + z + z',y)-p{t,x + z' 1 y)\ < \p(t,x + z + z',y)\+ \p(t,x + z',y)\ 

(A.10) < c(ps(t,x + z - y) +p s (t,x - y)). 

Observe that if \x — y\ < t 1 /", i.e. the diagonal regime holds for ps{t,x — y), we 
can then use the global upper-bound ( ps{t , x + z — y) +ps(t, x — y)) < tsjz, in the 
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control (1A.10I) so that 

D 12 (t,x,y) < cexp(ct u )p s (t,x 

< cexp(ct u )ps(t,x 


y) [ 


dz'\z'\ 

\z'\ d+a 


'\A>kt 1/o 



y)t 


-2+K 


If now I* — y\ > t 1 /®, we have for all given e £ (0,1), ps(t,x + z — y) < c e ps(t,x — y) 
if z B(x — y,e\x — y |). Indeed, \x — y + z\ > \\x — y\ — |^|| > e\x — y\. On 
the other hand, if z £ B(x — y,e\x — y |), \z\ > (1 — e)\x — y\ and \z\~ l ' d+a ' > < 
((1 — e)\x — 2/|) - ^ d+ “) < c E t~ 1 ps(t,x — y), up to a possible modification of c E . Thus, 


D 12 (t,x,y) < cexp(ct“)|i 2+ ^p s {t,x - y) + 

. , f dz f \z'\dz' 

c E ps{t,x-y / , ld+a / | nd+a + 

JB(x-y),e\x-v\)cn{\z\>iti/°‘} l~r + “ J\z’\<\tV« |*T + “ 

1 f , s , f \ Z '\ dz ' 1 

ps{i - [1+ z ~ v ) dz L si ,«. frw } 


b'l<4 4 

< cexp(ct ul )ps(t,x — y)t~ 2+ * . 


This proves that: 

(A.11) Di(t,x,y) < cexp(cf“)ps(f, x — y)t~ 2+ «. 


Let us observe that we have some continuity for the singularity w.r.t. the stability 
index, w.r.t. to the diffusive case, as far as the small jumps are concerned (see 
equation (I A. 71) j. The large jumps deteriorate the singularity. Precisely, 


\D 2 (t,x,y)\ < 


dz' 


\ z '\ d+a 

dz 

I p(t, X + z’ + z,y) - p(t, x + z', y)\ . 


From the previous arguments, splitting again the small and large jumps in the z 
variable w.r.t. the characteristic time scale t 1 / 0 , we get: 



I Pit, x + z' + z,y) ~ p(t , x + z',y)\ 



< - exp (ct ul )ps{t, x + z' 


y)- 


This therefore gives following the dichotomy used for the term D r2 (t,x,y): 


Q 

\D 2 {t,x,y)\ < ^ exp {ct u )p s {t,x - y), 

which together with (IA.11D indeed gives the bound \D(t, x, y)\ < # exp {ct u )ps{t, x— 

y)- 

Let us emphasize from |KM10j that the expansion (lA.dl) also holds in the strictly 
stable case. The previous control then gives the result, similarly to the discussion in 
the previous paragraph if s £ [0,f,;/2]. For s £ [ti/2,ti] we take the spatial adjoint 
and the associated contribution can be analyzed similarly. Eventually, this analysis 
extends to the remainder terms in (1A.4I) . 
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Appendix B. Two-sided Heat Kernel Estimates 

We derive in this section two-sided estimates for the density of X z />. Precisely 
we have the following Theorem. 

Theorem B.l (Two-Sided Heat Kernel Bounds). 

- Diffusive Case: Under (Ac), there exists c := c(/3,(Ao))> 1 s.t. for a given 
T > 0 and for all (x,y) £ (R d ) 2 , one has: 

• For d = 1: 

C" 1 f f t>/3\ ( \ X ~V\ 2 \ , , ^ 1 

< exp(—cT p ) exp I -c y/3 I + exp(-cT) exp I -c j r/3 j I f 

< Vx a (x, y) < 

Z T 

^{exp(cr«)e X p(-c-'!^) + exp(crA))exp(-c-‘{l!9£} A )}. 

(B.l) 


• For d = 2 and x ^ y\ 


^-|exp(-cT /3 )exp ^ (I log(c 1/2 )|I| x _ y |< c -i/2 T /3/2 + l) 

+ exp(-cr)exp^-c| ^ T ^ } 2 ~^ | 
< Px a (x, y) < 

Z T 

^|exp(cT^)exp j (| log(c~ 1/2 )|g [j; _ af |< c i/ 2^/2 + l) 

+ exp(cTw))exp }• 

(B.2) 


For d > 3 and x ^ y: 


-i / exp(-cT^) 
T/3|a;-y|rf- 2 


f exp(cT 2 ) 

T0|z - y| d - 2 6XP 
(B.3) 


exp I — c 


• - y\' 2 \ exp(-cT ) 

rp/3 J + T0d/2 


■ exp - 




rp/3 

< Px a (x, y) < 


_ 1 \x — y\ 2 \ exp(cT T +^^) 




TPd/2 


exp —c 


„_if|a;- 2 /| 


f F-nr i 
It/ 3 / 


2 ^ 

N 2-/3 


- Strictly Stable Case: Under (As), there exists c := c(/3,(As))> 1 s.t. for d > 2 
and a G (0,2) or d = 1, a < 1, a given T > 0 and for all (cc, y) £ (R d ) 2 ,x ^ y, one 
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has: 


.-if exp(-cr^) ¥ 

\ T 0 \x-y\d-^-y\< Tfi/a + 


1 exp(— cTP) 


< Px a (x, y) < 


/ exp(cT^) 1 exp(cr i -“<"-«)\ 1 

C \TP\x - y\d-a\ x -v\< TP,0 ‘ + jy3 d/a ~ N d.+a I ’ & — Al. 


(i + teS)' 


Q! 


(B.4) 

If d = 1, a > 1 we get for all (x, y) € 
-if 1 exp(-cT^) 


r pfid/a / A \d+a 


(i + fs* 1 )' 


.1 exp(cT 1 -“ j o-/3)) 

<Px z ,(x,y) < C <j y ^ /a ^ d+a h" 


( 1 + 


T@/ a 


1 

a 


(B.5) 


Proof. We first mention that the arguments in the proof of Lemma 14.81 also give 
that, a lower bound homogeneous to (14.611 holds for p z p(T,u). Precisely: 


(B.6) 




< p z e ( T,u ). 


- Diffusive case. Let us concentrate on the lower bounds, since the upper bounds 
are already derived in Corollary 14.11 and on the two-sided bounds for the two 
dimensional case. Note that we obtain in equations (IB. 21) . (IB. 81) . for dimensions 
d > 2, the same singularities that we observe for the Poisson kernel of the associated 
dimension. 

We focus below on the lower bound for d > 3 and the two-sided bounds for d = 2. 
The other controls can be derived from arguments similar to those developed below. 
One of the key points for the proof is the lower bound for the density of the spatial 
motion. Namely a lower bound homogeneous to the one in Lemma |4. 101 holds for 
the density of the spatial motion. It is classical, from chaining arguments, see e.g. 
Chapter 7 in Bass [Bas97| . to derive from the paranretrix expansions presented 
above that there exists c := c((A^)) > 1 s.t. for all (x,y) € ( R d ) 2 : 

(B.7) p(t , x, y) > c" 1 eX ^ /2 rt ' > exp ^ ^ . 

Lower bound for d > 3. We now start from (12.141) . (IB. 61) and (IB. 71) to derive: 


Px z p (x, y) > Tjyf | exp(-cT /3 ) ^ 


o m 


du , lx — y| 2 . 
■exp (W-—) 


d/2 


P + OO 


+ 


du 


It? u 


d/2 


exp (—c 


\x~y | 2 


) exp(— cu) exp I —c 


11/B-/3) 


(B.8) 


:= (mi + TO 2 )(T /3 , x — y). 


To control mi(T^,x — y) we consider again the previous, diagonal/off-diagonal 
dichotomy: 
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For \x — y\/T^ / 2 < 1 (diagonal regime) write: 

c -1 exp(— cTP) du 


mi 


{T p ,x-y) > 


TP 


,d/2 


> 


J\x-y\*/2 u 

c _1 exp(— cTP) du c _1 exp (—cT^) 

J\x-y\*/2 \x - y\ d ~ TP\x - y \ d ~ 2 


tp 


For |a; — y\/T& / 2 > 1 (off-diagonal regime) write: 


m i 


{T 0 ,x-y) > 


c 1 exp(— cTP) f T du f\x — y\ 


TP\x — y\ d 2 Jtp /2 u V u1 / 2 


d -2 


exp —c 


lz- y| 2 


> 


c 1 exp(—cT^) 
- y\ 


d—2 


■ exp ( —c 


> 


c " 1 exp(-cT /3 ) 
T^|cc — y\ d ~ 2 

We have thus established that: 


exp —c 


k~Z/| 

T& 

-y? 

tp 


2 ' /■ T '’ du 


/ TP /2 


(B.9) 


mi (T 13 , x — y) > 


c 1 exp(—cT^) 


■ exp —c 


T/ 3 


T/5|x-y| d - 2 

On the other hand for the contribution m 2 (TP,x — y) in (IB. 81) we get: 
For \x — y\/TPt 2 < 1 (diagonal regime) write: 


( T s w c 1 exp(-cT^) r 2T ‘ 3 
m 2 (T p ,a;-y)> —— / duexp ( -c 


u "i !/(!—z 3 )^ c 3 exp(— cTP) 


tp_ 


TP d /2 


T P(l+d/2) J T/} 

For \x — y|/T ^/ 2 > 1 (off-diagonal regime), write first from Young’s inequality, 

aScsK^w-w + frOjW]. 

Computations similar to those in Lemma 14.111 then yield: 


rmB \^ c exp(-cT) 
m 2 {T p ,x-y) > - —p - 


_ 3 1-3 

r 2[|z-y| 2 TT=7?]3=7? ^ 


l[\x-y\ 2 TT^P]?=Tt U d / 2 


exp —c 


\x - y| 2 


exp —c 


1/(1—Z 3 ) 


c 1 exp(—cT) / 

^ ^-exp -c 


\x-y \ 2 

tp 


2-/3 


-T/ 3 . 

1 


c 1 exp(— cT) ( 

> -= 7 ;-exp — c 


TP 


■~y ? 

tp 


\ X - y\ 2 TT&r ]*=f (rf/2-1) 


rb)l3^(«*/a-i) 


c 1 exp(—cT) / 

^-^75- exp ~ c 


[J>£(l+T=^)j 5=f 


|x-y | 5 


7 1 / 3 


2-/3 
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recalling that \/z > 1, exp(— z 1 ^ 2 ^)z 2 =%( d / 2 ^ > c 1 exp(— cz 1 ^ 2 for the 
last but one inequality. We have thus proved in all cases: 


m 2 (T /3 , x -y)> 


c 1 exp (—cT) 

rW 2 


■ exp —c 


\x-y I 2 

tp 


2=13 


which together with mm and (IB. 81) completes the proof of the lower bound for 
d > 3. 

Bounds for d = 2. The lower bound (IB. 81) and an homogeneous upper bound 
still hold, with obvious modifications of the constants. We focus on the contribu¬ 
tion toi(T /3 ,x — y) yielding the additional spatial diagonal singularity. The term 
to 2 (T^, 2 ; — y) can be analyzed as above. Let us write: 

c -1 exp(— cTP) rTlS 


TOi(T^, x — y) > 


TP 


du 

— exp [ —c 




Assume first that cf! 2 \x — y\/TP/ 2 < 1. Setting v := exp ^ in the above 

integral we derive: 


m\{T^, x — y) > 


c _1 exp(—cT^) dv 


TP 


- log(v) 


> 


c 1 exp(— cTP) 

T 


v log( log(u)) |q X p( Tli 


r exp(-e |a, T g l ) 


log(— log (v))dv 


> 


> 


c 1 exp(— cTP) 

Tp 


c 1 exp (—cTP) 


{ exp(- 


\x~y | 2 w , , \x-y \ 2 


-)(-l°g(c 


/•exp 

L 


TP n TP ^ 

l°g(— log(v))du| 


t p 


\x-y \ 2 


\x-y I 5 


tp 


{ exp(-c K T/ /' )( - log(c^)) 


Tp 


+ 1 og (c ^ X )(exp(-c |a: T/ f l ) — e a )}. 


-y | 2 


Hence: 

mi(T 0 ,x-y) > 

> 


c~ 1 e~ 1 exp(-cTP) exp 


TP 

The upper-bound could be derived similarly. 
If now cfl 2 \x — y\/T 0/2 > 1, we write: 


log(< 


\x-y | 2 

TP 


)l- 


mp 


(T^,x-y) > 


c 1 exp(-cTP) exp(~c ^ x T P ) /' T du 


TP 


I TP/ 2 u 


> 


c 1 exp(—cT^) exp(—c ^ T J^ ) 




The previous bounds give the result. 
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- Strictly Stable Case. We focus here on the lower bound for d > 2. The other 
cases can be handled similarly. From the lower bound in Kolkoltsov [KolOO] for the 
density and a chaining argument one derives that there exists c := c((As))> 1 s.t. 
for all t > 0, (x, y ) G (R d ) 2 : 


(B.10) 


p{t,x,y ) > 


c 1 exp(— ct) 


t <l/c 


1 + 


t i/° j 


d-\-a ' 


From (12.141) . (IB. 61) and (IB.101) we now get: 


c~ 1 r Tl3 

px z ^ x,y ) - yp \ ex p(~ cT/3 ) J o 


du 


,d/a 


r +oo 


du 


/ u 


d/ot 


(i + fe?) 


d+a 


exp(— cu) exp ( —c 


(i + fes 1 ) 

U 1 


d+a 


(B.ll) 


tp. 

:= (mi + m 2 )(T p ,x- y). 


Let us first control mi(T^,x — y) exploiting again the diagonal/off-diagonal di¬ 
chotomy. 


- If |a: — y\ < then 

,m0 w c” 1 exp(-c7+ [\ x ~ y \ a du c" 1 exp(-cT /3 ) 

mi{T - x - y) - — Tf — ^ ' 


- If |a; — y| > T^l a we get: 

C -1 rT p 

mi(T /3 ,x-y)>—g-exp(-cT f! ) du 

Jtp/2 


> 


c -i T P 


Tf >/2 \x-y\ d+a ~ \x-y\ d + a 


■exp (-cT p ). 


We have thus established 

(B ' 12) 1 T 13 

mi(T^x-y) > c- 1 exp(-c?+{ ^ _ y ^ d _J \ x -y\<T^+ |g _ y | d+a V»l>^/° }• 

Let us now turn to m 2 (T^, x — y). We get: 

- If \x - y\ < T^ a , 

r - 1 , r2T^ , -l 

m 2 (T^x-y) > _exp(-cT /3 ) J^ exp(-cT+ 

- If \x - y\ > T^/ a , 

c -i r' 2Td c~ 1 T ld 

m 2 (T^ ,x — y)> -yy— exp(— cT^) / udu > - -ry— exp(— cT 13 ). 

v - TP\x-y\ d+a v 1 Jtp - \x-y\ d+a 

Plugging the above estimates and (IB.121) into (IB. Ill) gives the result. 


□ 
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Remark B.l. Let us emphasize that the bounds of Theorem EH would hold under 
the weaker assumptions that the coefficients b,a are measurable and s.t. era* is 
Holder continuous and b is bounded. Indeed, the two-sided heat kernel bounds for 
the spatial motion hold in that case. We can for instance refer to Sheu |She91j in 
the diffusive case or to Huang lHua!5l in the strictly stable case. 

Let us point out as well that the two sided bounds hold in the diffusive case 
under those assumptions for the density of associated with the Euler scheme 

approximation for the spatial motion. Again, the key estimate is a two-sided bound 
for the Euler scheme which can be found in Lemaire and Menozzi |LM10j . In the 
strictly stable case, the upper bound holds under the indicated assumptions. This 
is a consequence of the parametrix expansion for the density of the scheme, see 
IKMIOj . The lower bound is more delicate to obtain since even in the diffusive 
case, the localization arguments which are standard for the SDE need to be carefully 
adapted for the scheme. 

Remark B.2. We conclude saying that the results in Theorem 13. II could be slightly 
improved in light of the sharp estimates of Theorem IB. II for d = 2. We did not 
exploit those controls in the presentation of the main results mainly for notational 
coherence and simplicity. 
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